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Abstract. A deep, fuly sampled diffraction limited (FWHM ~ 70 mas) narrow-band 
image of the central region in M87 was obtained with the Wide Filed and Planetary 
Camera 2 of the Hubble Space Telescope using the dithering technique. The Ha-|-[N ll] 
continuum subtracted image reveals a wealth of details in the gaseous disk structure 
described earlier by Ford et al. (1994). The disk morphology is dominated by a well 
defined three-arm spiral pattern. In addition, the major spiral arms contain a large 
number of small "arclets" covering a range of sizes (0"l-0"3 = 10-30 pc). The overall 
surface brightness profile inside a radius ~ l'.'5 (100 pc) is well represented by a power- 
law /(/i) ~ but when the central ~ 40 pc are excluded it can be equally well 
fit by an exponential disk. The major axis position angle remains constant at about 
PAdisk ~ 6° for the innermost ~ 1", implying the disk is oriented nearly perpendicular 
to the synchrotron jet (PAjct 291°). At larger radial distances the isophotes twist, 
refiecting the gas distribution in the filaments connecting to the disk outskirts. The 
ellipticity within the same radial range is e = 0.2 — 0.4, which implies an inclination 
angle of i ~ 35°. The sense of rotation combined with the dust obscuration pattern 
indicate that the spiral arms are trailing. 



1 Introduction 

The disk of ionized gas in the nucleus of M87 is currently the best example of a 
family of similar small (r ^ 100 pc) gaseous disks found to be common in the 
centers of elliptical galaxies with active nuclei (for a review see Ford et al. 1998). 
Several HST kinematical studies have shown that in M87 the gas is in Keplerian 
rotation, orbiting a massive black hole with a mass A/bh ^ 2 — 3 x IO^Mq 
(Harms et al. 1994; hereafter H94, Ford et al. 1996a, b; and Macchetto et al. 
1997, hereafter M97). The few other galaxies studied kinematically so far (NGC 
4261 - Ferrarese et al. 1996, NGC 6521 - Ferrrarese et al. 1998, NGC 4374 
- Bower et al. 1998) have further shown that nuclear gaseous disks offer an 
excellent tool for measuring the central black hole mass. 

Recent studies have revealed other important characteristics of the nuclear 
disk in M87. Its aparent minor axis (F96, M97) is closely aligned with the syn- 
chrotron jet {A9 ~ 10° — 15°) suggesting a causal relationship between the disk 
and the jet. The system of filaments in the center of M87 (Sparks, Ford & Kin- 
ney 1993; SFK) may also be causally connected to the disk. For example, the 
filaments extending ~17" (1200 pc) to the NW at PA ~ 315° are blue shifted 
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with respect to systemic velocity and show dust absorption implying they are 
on the near side of M87 as is the jet. These two findings led SFK to conclude 
that these filaments are streamers of gas flowing away from the center of M87 
rather then falling into it. The images in F94 (sec also Ford & Tsvetanov, this 
volume, FT98) show an apparent connection between at least some of the larger 
scale fiUaments and the ionized nuclear disk. 

Direct spectroscopic evidence for an outfiow was found recently. Several UV 
and optical absorption lines from neutral and very mildly ionized gas were mea- 
sured in the FOS spectrum of the nucleus (Tsvetanov et al. 1998; T98). These 
lines are broad (FWHM ~ 400 km s~^) and blue shifted by ~ 150 km s~^ 
with respect to M87's systemic velocity implying both an outflow and turbu- 
lence. In addition, non-circular velocity components - both blue and red shifted 
- were found at several locations in the disk (F96, FT98), and observed emis- 
sion lines are much broader than the expected broadening due to the Keplerian 
motion accross the FOS aperture. All these properties are best understood if a 
bi-directional wind from the disk were present. This wind may be an important 
mechanism for removing angular momentum from the disk to allow accretion 
through the disk onto the central black hole. 

Whatever the physical conditions in the disk it is important to map its mor- 
phology in detail. The first HST images (F94) have hinted that a spiral pattern 
could be present, but the signal-to-noise was too low for a definitive conclusion. 
In this paper we present deep, fully sampled diffraction limited narrow band 
images of the nuclear region in M87. We use these images to characterize the 
ellipticity, brightness distribution, and morphology of the disk. In this paper we 
adopt a distance to M87 of 15 Mpc, corresponding to a scale of 1" = 73 pc. 

2 Observations 

M87 was imaged though the WFPC2 narrow-band filter F658N (central wave- 
length/effective width = 6590/28 A) which isolates the Ha-|-[N ii] emission line 
complex at the systemic redshift. The nucleus was placed near the center of the 
planetary camera (PC) and the HST roll angle was 307°, orienting the M87 
synchrotron jet roughly along the diagonal of the PC field of view. 

The PC has a scale of 45.54 mas pixel"^ (Holtzman et al. 1995) and under- 
samples the HST diffraction limited PSF at all wavelengths shorter of ~ 1/i. 
To achieve full sampling we took images at four adjacent positions, offsetting 
the telescope by 5.5 PC pixels between each images - the so called dithering 
technique. The pattern used is illustrated in Fig. 1. At each of the four subpixel 
positions we took 3 images for cosmic ray rejection and cleaning of permanent 
defects. We re-calibrated the images using the best recommended calibration 
flies. 

To build a fully sampled image we used the variable-pixel linear reconstruc- 
tion algorithm (informally known as "drizzling") developed originally for the 
Hubble Deep Field (Fruchter & Hook 1997). The relative offsets of individual 
images were estimated from the positions of good signal-to-noise objects in the 
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Fig. 1. Left panel illustrates the dithering pattern. Three images were taken at each 
of the positions numbered 1 to 4. The telescope offset was 5.5 PCI pixels between 
positions. The right panel shows the radial profile of an unsaturated star ~6" south of 
nucleus. The overplotted Gaussian has a FWHM equal to the HST diffraction limit at 
the wavelength of observation (6580 A). The first Airy ring is marked by the arrow. 



field and also through a 2-D cross correlation. Both techniques yielded similar 
results. No significant offsets were found for the three frames taken at each one 
of the subpixcl positions. These were combined to remove the cosmic ray events. 
The four combined frames were then combined on an output grid of 25 mas 
pixel" ^ using the drizzling algorithm. The 25 mas pixel size is a good match 
for the ^ 70 mas HST diffraction limit at Ha. This is illustrated by the radial 
profile of a star in the field shown in Fig. 1 where the first Airy ring is clearly 
seen. 

Because of serious observing time restrictions, and also because the underly- 
ing galaxy profile is known to be smooth, no continuum images were obtained 
simultaneously with the on-line ones. To create a suitable off-band image for 
continuum subtraction we used a number of M87 images from the HST public 
archive. These include a series of F814W images taken on November 11, 1995, 
which have small rcaltivc offsets that allow us to build a fully sampled image 
using the same drizzling technique. The offset pattern in not ideal (as in the case 
of F658N), but this is less critical because of the smooth shape of the continuum 
light distribution. 

To form an image suitable for continuum subtraction at the wavelength of 
Ha we used a {V — I) color image to correct the F814W image for color effects, 
which are particularly important at the positions of the dust bands. Finally, the 
continuum subtracted Ha-|-[N ii] image was flux calibrated using the estimated 
system throughput at the redshifted position of Ha and the averaged mesured 
[N ii]/Ha line ratio. 
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3 Disk morphology 

The central 1" — 2" of the continuuin subtracted Ha+[N ii] image (Fig. 2, left 

panel) is dominated by a clockwise winding 3-arm spiral pattern superposed 
on an underlying disk-like morphology. At larger radial distances the gaseous 
structure is less well organized, asymmetric, and gradually transition into the 
larger scale Ha filaments. Several of these filaments can be traced to connect to 
the disk but none appears to go directly into the nucleus. 




-a -1 1 -2-10 1 

offset, arcsec Offset, arcsec 



Fig. 2. The Hq-|-[N ii] continuum subtracted image is shown in the left panel. The 
gray scale saturates the central few pixels, which are uncertain due to the imperfect 
color matching and variability of the central point source. The right panel shows the 
ratio of the observed structure to the smooth disk model, as duscussed in the text. The 
image is stretched linearly between ±70% relative to the model. 

To estimate the parameters describing the surface brightness profile, elliptic- 
ity and orientation, we have performed an elliptical isophote analysis. Wc started 
by setting all parameters free, and gradually introduced constraints such as fixed 
center position, major axis position angle and/or ellipticity, in an attempt to ex- 
plore the robustness of the solution. To our satisfaction the values for all major 
parameters of the fit remained stable regardless of the constraints used. 

Fig. 3 shows the surface brightness radial profile, major axis position angle, 
and ellipticity as a function of radial distance. Outside r ^ 0"5 (40 pc), the sur- 
face brightness profile is well described by an exponential disk with an effective 
radius of 45 pc. There is, however, a significant excess of light in the central 0'.'5, 
and a single power law /pi ~ (r/50)~^-^^ is a better representation of the overall 
profile. In the region 20-90 pc, the orientation of the major axis remains roughly 
constant at PAmaj '--^ 6° . At the same time the ellipticity increases smoothly from 
0.2 to 0.4. The constancy of PAmaj is usually interpreted as a signature of a disk, 
while the ellipticity could be infiuenced by the spiral arms or there could be a 
warp. The central few points are less reliable because of possible color mismatch 
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Fig. 3. Left: The surface brightness profile of the M87 nuclear gaseous disk. The short 
dashed line is a power law 7pi ~ (r/50)~^'^^, and the long dahsed line is an exponential 
disk with /disk ~ exp{—r/45), where r is in pc. Right: The radial dependence of the 
major axis position angle (PAmaj) and ellipticity (e). 

of the on- and ofF-band images and the variability of the nuclear point source 
(T98). 

The kinematic studies (F96, M97) have shown that the north part of the disk 
is receiding and the south part is approaching and the estimated inclination is 
consistent with the one inferred from isophote analysis. The obscuration patches 
on the SE side (F94) indicate that the SE side is closer to us. Combining all this 
information indicates that the spiral arms are trailing. 

This work is supported by NASA grant NAG-1640 to the HST FOS team. 
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